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ABSTRACT 

We present a study of MSX point sources in the Galactic bulge < 3°, 1° < \b\ 
< 5°), observed at A. C, D and E-band (8 to 21 /im), with a total area ~ 48 deg 2 and 
more than 7000 detected sources in the MSX D-band (15 /jm). We discuss the nature 
of the MSX sources (mostly AGB stars), their luminosities, the interstellar extinction, 
the mass-loss rate distribution and the total mass-loss rate in the bulge. The mid- 
infrared data of MSX point sources have been combined with the near-infrared (J, 
H and K s ) data of 2MASS survey. The cross-identification was restricted to A' s -band 
detected sources with K s < 11 mag. However, for those bright MSX D-band sources 
([D] < 4.0 mag), which do not satisfy this criteria, we have set no AVband magnitude 
cut off. The bolomctric magnitudes and the corresponding luminosities of the MSX 
sources were derived by fitting blackbody curves. The relation between M and (K s - 
[15])o was used to derive the mass-loss rate of each MSX source in the bulge fields. 
Except for very few post-AGB stars, planetary nebulae and OH/IR stars, a large 
fraction of the detected sources at 15 /jm (MSX D-band) are AGB stars well above 
the RGB tip. A number of them show an excess in ([A]-[D]) and (A S -[D]) colours, 
characteristic of mass- loss. These colours, especially (A s -[D])o, enable estimation of 
the mass-loss rates (M) of the sources in the bulge fields which range from 10 to 
10~ 4 M yr _1 . Taking into consideration the completeness of the mass-loss rate bins, 
we find that the contribution to the integrated mass-loss is probably dominated by 
mass-loss rates larger than 3 x 10~ 7 M yr _1 and is about 1.96 x 10~ 4 M Q yr _1 deg -2 
in the "intermediate" and "outer" bulge fields of sources with mass-loss rates, M > 
3 x 10 _7 M Q yr _1 . The corresponding integrated mass-loss rate per unit stellar mass 
is 0.48 x 10 -11 yr _1 . Apart from this, the various mid- and near-infrared colour-colour 
and colour-magnitude diagrams are discussed in the paper to study the nature of the 
stellar population in the MSX bulge fields. 

Key words: stars: AGB and post-AGB - stars: circumstellar dust - stars: mass-loss 
- dust: extinction - infrared: stars 



1 INTRODUCTION 

Study of the stellar population of the Galactic bulge fields 
is of prime importance and plays a crucial role in under- 
standing the formation and evolution of the Galaxy. The 
high luminosity Asymptotic Giant Branch (AGB) stars are 
ideal tracers of stellar population in regions of high ex- 
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tinction such as the "intermediate" and "outer" Galactic 
bulge. Intense mass-loss ( > 10 -6 M© yr _1 ) phase has been 
identified with stars evolving along the AGB phase. Hence, 
these sources are enshrouded with circumstellar envelopes of 
dust and gas. The low effective temperatures and thermal 
emission from warm dust make these stars bright in the in- 
frared. Deep and large area infrared surveys offer a unique 
view of the stellar population towards the inner Galaxy 
as high interstellar extinction hinders the study at opti- 
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cal wavelengths. A combination of mid- and near-infrared 
data is essential to sample the high mass-loss AGB popu- 
lation. Mid-infrared data are more sensitive to the infrared 
excess which is a consequence of mass-loss in the AGB stars. 
With the availability of infrared data from surveys like DE- 
NIS (Epchtein et al. 1994), 2MASS (Beichman et al. 1998; 
Skrutskie et al. 2006), ISOGAL (Omont et al. 2003) and 
MSX (Price et al. 2001), there have been a large number 
of studies on the AGB population of the Galactic bulge 
fields (e.g. Glass et al. 1999; Schultheis & Glass 2001; Glass 
& Schultheis 2002; Groenewegen & Blommaert 2005). In a 
recent study, Ojha et al. (2003) have combined the mid- 
and near-infrared photometry from ISOGAL, DENIS and 
2MASS to study the nature of ISOGAL sources in the "in- 
termediate" Galactic bulge and discuss their mass-loss rates. 
Further, the data from GLIMPSE II survey (Ed Churchwell; 
private communication) will provide an excellent opportu- 
nity to study the AGB stars in the inner bulge region. How- 
ever, the sample has to be restricted to high extinction re- 
gions owing to saturation effects of the survey. In the present 
study we have not included GLIMPSE II data as it does not 
cover the entire "intermediate" and "outer" bulge regions 
selected in this paper. The stellar population study of the 
inner bulge with GLIMPSE II data will be presented in a 
future paper. The AGB stars contribute more than 70% to- 
wards the enrichment of the dust component of the ISM in 
the solar neighbourhood (Sedlmayer 1994) and hence it is 
important to study their mass-loss in different parts of the 
Galaxy. As discussed in Ojha et al. (2003), the total mass 
returned to the ISM is dominated by mass-loss rates greater 
than 10 -6 Mq yr _1 and hence the detection of the entire 
population of the high mass-loss stars becomes important 
for determination of the total mass returned to the ISM and 
probably the total mass of the bulge. 

In this paper, we report the study of the MSX point 
sources in the "intermediate" and "outer" Galactic bulge 
fields (| I | < 3°, 1° < |6| < 5°) with a total area ~ 48 
deg 2 and more than 7000 detected sources in the MSX D- 
band (15 /mi). Here, the division into the two aforemen- 
tioned bulge fields is based on the Galactic latitude. The 
MSX bulge fields in the Galactic latitude bins 1° < |6| < 2° 
and 2° < |6| < 5° are defined as the "intermediate" and the 
"outer" bulge regions, respectively. We restrict our study 
to the "intermediate" and "outer" Galactic bulge because 
the extinction here is much less and more homogeneous as 
compared to the inner bulge. The density of the sources in 
"intermediate" and "outer" bulge also reduce the number of 
spurious associations to the minimum. As shown in Price et 
al. (2001), uniform sky coverage in the MSX bands extends 
for only the northern latitudes upto b — 6° and the higher 
latitudes (outer Galaxy) are very sparsely sampled. We have 
restricted the latitude selection such that the fields covered 
by us are free from the non-uniform sky coverage beyond 
|6| > 5°. We have combined the mid-infrared data of the 
MSX point sources with the near-infrared 2MASS data to 
determine their nature and the interstellar extinction. Most 
of the sources are AGB stars well above the RGB tip with 
high mass-loss. We have determined the luminosities and 



mass-loss rates of the stellar population in the MSX bulge 
fields. 

The outline of the paper is as follows : In §2, we present 
the MSX and 2MASS observations and describe the cross- 
identification between MSX and 2MASS sources. In §3, we 
discuss the determination of interstellar extinction in the 
line of sight of the bulge fields using the isochrone fitting 
method. The ISOGAL and DENIS associations of the MSX 
sources in the bulge fields are discussed in §4. §5 presents 
the derivation of bolometric magnitude (Mtoi) and luminos- 
ity for each star in the bulge fields. In §6, we discuss the 
derivation of mass-loss rates based on (if 8 -[D])o colour and 
the estimation of the total mass-loss rate in the "interme- 
diate" and "outer" bulge. In §7, we present the nature of 
the MSX sources from the various colour-colour and colour- 
magnitude diagrams. We summarize our results in §8. 



2 OBSERVATIONS AND 

CROSS-CORRELATION OF MSX AND 
2MASS SOURCES 

We have used the deep MSX Point Source Catalog Ver- 
sion 2.3 (Egan et al. 2003) in this paper. The catalogue 
(MSXPSC V2.3) lists the sources detected in MSX mid- 
infrared bands A, C, D and E with A(AA) corresponding 
to 8.28(3.36), 12.13(1.72), 14.65(2.23) and 21.34(6.24) /xm, 
respectively. MSXPSC V2.3 has several improvements over 
the initial published catalogue, MSXPSC Version 1.2 (Egan 
et al. 1999). In this latest version, the photometry is based 
on co-added image plates, as opposed to single-scan data, 
which results in improved sensitivity and hence reliability 
in the fluxes. Comparison with Tycho-2 positions indicates 
that the positional accuracy, a — 2" at 8 /im, of the new 
catalogue is better as compared to MSXPSC VI. 2. Also, 
MSXPSC V2.3 is 100% complete in the A-band (8 /um) and 
approaches 100% completeness in the other three bands at 
the survey sensitivity limit of the MSX image data (~ 0.1 
Jy in A-band, which is more sensitive than the other bands 
by a factor of ~ 10) (Egan et al. 2003). 

The total number of sources detected in our fields are 
29709, 8125, 7390, and 3380 in the MSX A-, C-, D- and 
E-bands, respectively. In the study that follows, we have in- 
cluded only good quality MSX data for which the flux qual- 
ity flags are 3 and 4 (which implies S/N > 7). Taking the 
quality flags into account, the total number of good quality 
sources in our fields are 24035, 2571, 2550 and 866 in the 
four bands, respectively. The corresponding average source 
densities are ~ 5.0xl0 2 deg" 2 , ~ 0.5xl0 2 deg~ 2 , ~ 0.5xl0 2 
deg -2 and ~ 0.2xl0 2 deg -2 in the four bands, respectively. 

The MSX catalogue is more sensitive (by a factor of ~ 
4 in the A-band) towards the inner Galactic longitudes as 
compared to the outer longitudes. This is due to the fact 
that additional maps of regions at I = were co-added to 
the survey data (Sean Carey; private communication) . These 
additional maps are ~ 4 times more sensitive. This is evident 
from the inspection of the data sets. Hence, for clarity we 
have divided the catalogue into two sets for our study: "in- 
ner" zone < 0.5°) and "outer" zone (0.5° < \l\ < 3.0°). 
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Figure 1. A, C, D and E-band source distributions in half magnitude bins for good quality (quality flags 3 Bz 4) MSX data. Magnitudes 
are converted into fluxes using the zero magnitude flux given by Egan et al. (1999). The solid line denotes the MSX sources for the 
"inner" zone < 0.5°) and the dotted line represents the sources for the "outer" zone (0.5° < \l\ < 3.0°) in the bulge. 



The magnitude histograms of the two zones in the MSX A-, 
C-, D- and E-bands are displayed in Fig. [TJ The histograms 
show the total source counts of the good quality sources. 
The 50% completeness limits (as determined from the flux 
histograms) of the "inner" zone are 0.08, 0.40, 0.32 and 1.00 
Jy for A-, C-, D- and E-bands, respectively. For the "outer" 
zone, the 50% completeness limits correspond to 0.16, 1.58, 
1.26 and 3.98 Jy for A-, C-, D- and E-bands, respectively. 

From the ISOGAL survey, it is seen that a number of 
the ISOGAL 15 /im sources are AGB stars in the Galac- 



tic bulge and the central disk with evidence of mass-loss in 
majority of the cases (Omont et al. 2003). The mass- loss 
in AGB stars may be characterized by their 15 (im excess 
(Omont et al. 1999; Ortiz et al. 2002; Ojha et al. 2003). The 
MSX D-band is close to this wavelength and hence in the 
study that follows, our prime focus lies on the sources de- 
tected in the D-band. The proportions of D-band sources 
associated with A, C, and E-band sources are ~ 100%, ~ 
91% and ~ 33%, respectively. 

The near-infrared data used in this paper have 
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been obtained from the 2MASS Point Source Catalogue 
(2MASSPSC) in the three bands, J (1.25 (im), H (1.65 fim) 
and K s (2.17 /Ltm). The full sky 2MASS catalogue covers the 
entire MSX bulge fields. The 2MASS PSC is >99 % com- 
plete in the absence of confusion at the IOct sensitivity limit 
of 15.8, 15.1 and 14.3 mag in the J-, H- and A s -band, re- 
spectively (Skrutskie et al. 2006). However, the turnover in 
the source counts in the Galactic plane field occurs nearly 1 
mag brighter, because of the effects of confusion noise on the 
detection thresholds. The primary areas of confusion are (1) 
longitudes ±75 deg from the Galactic center and latitudes 
1 deg from the Galactic plane and (2) within an approxi- 
mately 5 deg radius of the Galactic center (Skrutskie et al. 
2006). 

A cross-correlation algorithm similar to the one used by 
Schuller et al. (2003), to associate the ISOGAL 7 and 15 
sources with DENIS sources, has been used to search for 
2MASS counterparts of the MSX point sources in the bulge 
fields. However, to restrict the number of spurious detections 
at the fainter end, we set our magnitude cut at conservative 
brighter level and retain only those 2MASS sources with 
K$ < 11 mag for cross-correlating with the MSX catalogue. 
This 2MASS magnitude cut (K a < 11 mag) corresponds to 
an average density of A^-band sources n ~ 22 300 per deg 2 . 
However, the density of K B < 11 mag sources varies from 
18 030 to 26 935 per deg 2 depending on the locations of 
the bulge fields. To search for 2MASS counterparts, we as- 
sumed a main association radius (r a ) of 4.0". The density 
limit and main association radius are chosen such as to re- 
strict the chance association to less than 10% (where chance 
association, y = n7rr 2 < 0.1). As a result, ~ 89% of the MSX 
D-band sources (2265 sources with quality flags 3 and 4) 
within the 2MASS observations have an association with a 
2MASS source (K s < 11.0 mag). Also, 2% (48 sources) of 
the MSX D-band sources have a 2MASS counterpart which 
are saturated (K a < 3.5 mag). In Fig. [5] we show the 3-D 
plots of the MSX source density, the mean Av and the total 
number of 2MASS sources (K s < 11.0 mag) in the various 
subfields of the bulge as a function of the Galactic longitude 
and latitude. This information is also presented in tabular 
form in Table fl 

Keeping in mind the importance of high mass-loss AGB 
stars in our study, we introduce two additional samples here. 
We search for 2MASS counterparts of bright MSX sources 
detected in the D-band but not associated with 2MASS 
sources with A' s < 11.0 mag. We restrict our sample to 
sources brighter than 4th magnitude in the MSX D-band. 
This magnitude limit was chosen such as to include all 
the sources with high mass-loss rates (logio(M)(M0 yr _1 ) 
> — 6.0; see §7). There are 165 such sources with good 
quality MSX D-band photometry. For identifying 2MASS 
counterparts in this additional sample, we extend the main 
association radius to a conservative value of ~ 5" without 
any A s -band magnitude cut. The increase in the association 
radius by 1" enables us to include some known high mass- 



1 Table 1 is only available in the electronic form via the VizicR 
Service at the CDS. 



loss peculiar stars (e.g. carbon stars) in the sample (see §7). 
However, a crucial point is that the removal of the K s < 
11.0 magnitude cut increases the 2MASS source density to 
~ 10 5 sources per deg 2 . Hence, for a main association radius 
of 5", the chance association becomes very large (y « 0.6) 
compared to the 10% limit set for the rest of the sources 
in the bulge fields. The corresponding pure Poisson chance 
association is given by 1 — exp(— y) ~ 0.45. This implies 
that about half of the 2MASS associations are likely to be 
spurious. However, it should be kept in mind that all such 
spurious associations would possibly have a A a -band coun- 
terpart fainter than the 2MASS limit and hence in our study 
these sources are retained as a lower limit for A s -band and 
limit for the high mass- loss end (see §6). We are now left 
with two samples of these bright MSX D-band sources - 
1) Sources having a 2MASS counterpart within 5" radius, 
which we designate as sample-A and 2) Sources having no 
2MASS counterpart within this radius, which we name as 
sample-U. In sample-A, we have 125 sources out of which ~ 
70% have a single 2MASS association within the 5" radius. 
Sample-U contains 42 sources, and we assign a value of K B 
= 13.7 (which is the A s -band limit of the sample-A sources) 
for all these sources. 

For this study, alongwith the good quality MSX band 
data, we have used only good quality 2MASS J-, H- and 
A s -band magnitudes with 'rd-fiag' values between 1 and 3 
which generally implies best quality detection, photometry 
and astrometry. However, the 2MASS sources with 'rd-flag' 
value of 6 in the A s -band (which corresponds to a positive 
detection with an upper magnitude limit) were also included 
with good quality sources ('rd-flag' 1-3) for the determi- 
nation of mass-loss rate only (see §6). 

In Table we present the full catalogue of MSX- 
2MASS sources from the bulge fields (see Fig. [2] and Ta- 
ble 1), a sample of which is given in this paper. Figure [3] 
shows the histogram of positional difference of MSX D- 
band-2MASS cross-identified sources in the bulge fields. The 
mis of the differences of MSX-2MASS association is ~ 0.7". 
The histogram peaks at ~1.0" and the almost exponen- 
tial decreasing trend with increasing distance indicates good 
quality associations, which implies that the number of spu- 
rious associations would be minimal in our cross-correlated 
sample. 



3 INTERSTELLAR EXTINCTION AND 
FOREGROUND DISK STARS 

In most parts of the Galactic bulge, the interstellar ex- 
tinction is not homogeneous and occurs in clumps, hence, 
a detailed extinction map is essential in stellar population 
studies (Schultheis et al. 1999). Inspite of the recent im- 
provements in the extinction measurements of the Galactic 
bulge (Schultheis et al. 1999; Dutra et al. 2001; Marshall 
et al. 2006; Indebetouw et al. 2005), the determinations are 
still uncertain. In this paper, the method as described in 

2 Table 2 is available in the electronic form via the VizieR Service 
at the CDS. 
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Figure 2. The figure shows the 3-D plots of the MSX source density, the mean Ay and the total number of 2MASS sources (K B < 
11.0 mag) in the various subfields of the bulge as a function of the Galactic longitude and latitude. This information in tabular form is 
presented in Table 1 which is only available in the electronic form via the VizicR Service at the CDS. 



Table 2. Sample catalogue of MSX-2MASS sources in the bulge fields. 



Seq. 


Name (MSX6C-) 


RA (2000) 


I 


J 


H 


K a * 


[A] 


[C] 


[D] 


[E] 


A v 








Dec (2000) 


b 






K s ** 


[7] 




[15] 












(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


6007 


G359.9129- 


-02.1610 


264.2674 


-0.0871 


10.11 


8.16 


6.99 


4.15 


3.13 


2.96 


1.99 


9.48 








-27.8677 


2.1610 






7.11 


4.72 




99.99 






6057 


G359.9338- 


-02.1591 


264.2818 


-0.0662 


7.35 


5.81 


5.08 


3.94 


3.04 


2.91 


2.24 


3.98 








-27.8510 


2.1591 






6.31 


4.34 




99.99 






6444 


G359.9392- 


-02.0232 


264.4150 


-0.0608 


6.53 


5.08 


4.30 


3.45 


2.57 


2.51 


2.12 


3.53 








-27.9193 


2.0232 






99.99 


3.39 




99.99 






6838 


G000.3366- 


-02.1338 


264.5491 


0.3366 


9.46 


7.94 


7.03 


4.49 


3.50 


3.33 


2.50 


5.50 








-27.5248 


2.1338 






6.93 


99.99 




3.74 






6895 


G000.1512- 


-01.9967 


264.5683 


0.1512 


12.05 


9.72 


8.12 


3.95 


2.75 


2.56 


1.70 


14.59 








-27.7546 


1.9967 






8.33 


99.99 




2.45 







* 2MASS K B magnitude; ** DENIS K B magnitude; Magnitude assigned to value of 99.99 means that the source is not detected 
in the corresponding magnitude band. 



Ojha et al. (2003) which is based on the procedure outlined 
in Schultheis et al. (1999) is used to determine the inter- 
stellar extinction. The MSX catalogue was divided into 60 
subfields. This enables in decreasing the effects of variable 
and patchy extinction present in the entire bulge field. For 
\l\ < 1.0°, we have 0.5 deg 2 fields with steps of 0.5 deg in 
longitude and 1.0 deg in latitude. For rest of the fields, the 
field size was increased to 1.0 deg 2 with steps of 1.0 deg both 
in longitude and latitude. This ensures appreciable statisti- 



cal sample in each field to derive mean Ay. The fact that 
we use 1 deg 2 subfields to derive extinction gives some dis- 
persion of the order of typically 1.5 mag in Ay (see Table 
1). However, at Galactic latitudes greater than 1 deg (where 
our fields are located), the dumpiness of interstellar extinc- 
tion is less prominent (Schultheis et al. 1999). The subfields 
are listed in Table 1 and the field centres are displayed in 
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Figure 3. Histogram of positional difference (in arcsec) between 
good quality MSX D-band (quality flags 3 and 4) and 2MASS 
cross-identified sources. The MSX— 2MASS association was lim- 
ited to r < 4" for bulk of the sources and increased to 5" for 
sample- A sources (see the text). 
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Figure 4. Color-magnitude diagrams (J — K B /K B ) of 2MASS 
sources in a few selected MSX bulge fields. The CM diagrams 
for all the bulge fields will be available in electronic form via the 
VizieR Service at the CDS. 



the various plots of Fig. [2] Figure [jjj shows the J — K s /K s 
colour-magnitude diagrams for a few selected MSX fields 
presented in Fig. [2] and listed in Table 1. Figure [5] shows the 
J — K s /K s colour-magnitude diagrams for two sample fields 
(-0.5° < I < 0.0°, 3° < b < 4° ; -2.0° < I < -1.0°, -2° < 
b < -1° ). 

As is clearly seen in Figs. [4]and[5] the J—K s /K s colour- 
magnitude diagrams of 2MASS sources in the bulge fields 
show a well-defined red giant and AGB sequence shifted 
by fairly uniform extinction, with respect to the reference 
K a0 vs. (J - K a ) of Bertelli et al. (1994) with Z = 0.02 
and a distance modulus of 14.5 (distance to the Galac- 
tic Centre : 8 kpc). However, a point to be noted here 
is that the intrinsic depth of the bulge is about 0.3 mag 
(Glass et al. 1995). We have assumed Aj/Av = 0.256 and 
Ak s /Av = 0.089 (Glass 1999) and calculated the individ- 
ual extinction values of the 2MASS bulge sources as de- 
scribed in Ojha et al. (2003). The mean Ay value for each 
field has been determined by a Gaussian fit to the Ay 
distribution. It should be noted here that owing to low 
star counts, we have increased the size of some bins (e.g 
0.0 < I < 0.5; 4.0 < b < 5.0 & 0.5 < I < 1.0; 4.0 < b < 5.0 
are merged into one bin) to estimate the mean extinction 
values. The details of the subfields, giving the mean Ay, the 
total number of sources and the statistics of the background 
and foreground population are listed in Table 1. To derive 
the dereddened magnitudes for the MSX sources, we have 
used the extinction law - A A /A V = 0.022; A c /A v = 0.023; 
Ad/ Ay = 0.013; A E /A V = 0.016 (Messineo et al. 2002). 

As seen in Fig. [5] the 2MASS sources with anomalously 
low values of Ay are probably foreground stars. For each 
field we empirically define an isochrone "F" for which we 
assume that all the sources left of it are foreground stars. 
Also seen in Fig. [5] and in each of the bulge fields, these 
foreground stars are around the isochrone with Ay ~ mag 
and clearly left of the bulk of the stars grouped around the 
isochrone with mean Ay of the field. These foreground stars 
will be no longer considered in the following discussions of 
bulge stars. However, this empirical way of rejecting fore- 
ground population does not exclude foreground stars with 
significant mass- loss. 

There are also a number of stars with J — K s values 
significantly larger than the bulk of the other stars in each 
bulge field (right of the isochrone empirically defined as "B" 
in Fig. ■ The sources to the right of this isochrone are 
termed as "B-sources" . We can see three reasons for such an 
excess in J — K s : 1) intrinsic (J — K s )o excess induced by 
a large mass-loss, which should be accompanied by a large 
15 /im excess; 2) spurious association or wrong photometry 
which is rather unlikely for 2MASS sources well above the 
detection limit. 3) excess in Ay which should probably be 
due to a patchy extinction on the bulge line of sight (addi- 
tional extinction from dust layers behind the Galactic centre 
for background stars appears unlikely at such high galactic 



3 The J — K s /K s colour-magnitude diagrams of 2MASS sources 
in the bulge fields are available in electronic form via the VizieR 
Service at the CDS. 
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latitudes) , however this does not seem a major source of ex- 
cess with relatively small average extinction. For the bright 
stars in this region, with K a < 8 mag, we use the mean ex- 
tinction value (typical SAv ~ 1.56 mag) of the field. This 
is because of two reasons: many of the stars with K B < 8 
mag should have a large mass-loss which could produce an 
excess in ( J — A" s )o; and at the bright end we are restricted 
by the limit of the Bertelli et al. (1994) isochrone which is 
used as the reference in shifting the sources and determin- 
ing the individual extinctions. For rest of the sources in this 
region which include the faint sources, with K a > 8 mag, 
we determine their specific extinction from the J — K s /K s 
colour-magnitude diagram and the zero-extinction isochrone 
similar to the extinction determination for the bulk of the 
bulge sources. For few bright sources in the region left of 
isochrone "B", we have extrapolated the Bertelli isochrone 
to estimate the individual extinction values. 

The uncertainties in the estimated values of Av due 
to the assumption of a 10 Gyr stellar population of solar 
metallicity as the reference system has been discussed in 
Schultheis et al. (1998; 1999). It is shown that the age-range 
(of 5 Gyr) has hardly any effect, whereas, the assumption 
of solar metallicity results in typical uncertainties of ~ 1 
mag in the determination of Av- Another important point 
to be noted here is that for the high mass-loss AGB stars, 
post-AGB sources and planetary nebulae population, the 
contribution from the intrinsic extinction due to circumstel- 
lar dust has not been considered. This will result in larger 
uncertainties in the estimated Av values. 

The sample of MSX sources discussed in this paper 
contains ~ 24000 2MASS sources which are detected in 
the A s -band. From this sample we have identified "fore- 
ground" and "B-sources" based on the method described 
above and demonstrated in Fig. [5] Their fractions are about 
12% (2937) and 9% (2113), respectively. For the sample-A 
population, we have 111 MSX sources (after foreground sub- 
traction) which have good quality D magnitudes (D-band 
quality flags 3 & 4) and A s -band photometry ('rd-flag' val- 
ues 1 - 3). It should be kept in mind that for sample-U 
sources it is not possible to identify the foreground popula- 
tion and hence, they are not removed from the sample. 



4 ISOGAL AND DENIS ASSOCIATION 

We have also cross-correlated all the MSX sources with the 
ISOGAL sources in the bulge fields. The ISOGAL Point 
Sources (IGPSC; Schuller et al. 2003) have been obtained 
from the VizieR Service at CDS. ISOGAL Point Source Cat- 
alogue contains about 10 s stars, with associated K, J, and 
I DENIS data for most of them (see Schuller et al. 2003). 
The MSX and ISOGAL associations were searched within 
a radius of 6" (size of the ISO pixel). There are a total of 
about 2110 ISOGAL sources at 7 and 15 /im in the MSX 
bulge fields (total area ~ 48 deg 2 ). Out of this 194 ISO- 
GAL sources at 7 /im have a MSX A-band association and 
55 ISOGAL sources at 15 fim have MSX D-band associa- 
tion. The corresponding magnitude limits for these associ- 
ated sources are 7.8 and 4.6 mags in the A- and D-band, 
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Figure 5. Colour-magnitude diagrams (J — K B /K B ) of MSX- 
2MASS sources in two MSX bulge fields. The four isochroncs 
(Bertelli et al. 1994), placed at 8 kpc distance for a 10 Gyr pop- 
ulation with Z=0.02 are shown in the figure for Ay = mag, 
for the Ay limit adopted for the "foreground" sources (defined as 
"F"), for mean Ay value of the field and for the Ay limit adopted 
for the "B-sourccs" , respectively. 

respectively. We have checked the consistency of MSX and 
ISOGAL 15 /jm magnitudes for the 38 strongest sources 
(with [D] < 4.0 mag; among our sample-A). The average 
difference is [15]msx - [15]/so = -0.10±0.43 mag, where the 
large dispersion is probably related to the source variabil- 
ity. We have searched the DENIS counterparts for only the 
ISOGAL cross-correlated MSX sources. More than 99% of 
the ISOGAL sources in the bulge fields have a DENIS coun- 
terpart. There are 435, 61 and 54 2MASS A s -band sources 
which have associations with DENIS K s sources, DENIS A s 
& MSX D-band sources and DENIS K a , MSX D-band & 
ISOGAL 15 /im sources, respectively. For the ISOGAL cross- 
correlated sources, we have complemented the 2MASS and 
MSX sources in our catalogue (Table 2) by the DENIS and 
ISOGAL data. Within the ISOGAL fields (Ojha et al. 2003), 
more than 95% of MSX sources having 2MASS association 
have an ISOGAL counterpart. We have used the DENIS and 
ISOGAL data, wherever available, in the estimates of the 
bolometric magnitudes (see §5) and the mass-loss rates (see 
§6) in order to mitigate the effects of variability. Hence, for 
the bulge sources with DENIS and ISOGAL associations, we 
have used the average of the 2MASS and DENIS values for 
the As-band magnitude and the MSX and ISOGAL 15 /jm 
values for the D-band magnitudes. 



5 BOLOMETRIC MAGNITUDES AND 
LUMINOSITIES 

After dereddening (see §3), the bolometric magnitudes 
(Mboi) are derived by integrating the flux densities over the 
wavelength range between 1.25 /im < A < 15 /im by fitting a 
blackbody curve. For each object we used the dereddened J- 
, A s -, A-, D-band magnitudes and the estimated Av values 
as input parameters. The main error of the bolometric cor- 
rections results from the interstellar extinction values which 
gives an error of ~ 0.2 - 0.3 mag in Mbol- The assumption 
of a blackbody curve for the AGB spectra is less than opti- 
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mal given the presence of strong molecular bands. We have 
compared the results obtained from the blackbody fitting 
with that derived using multi-band bolometric correction 
for AGB stars (Cecile Loup; private communication). How- 
ever, this method is valid for limited range of colours and 
luminosities. The difference in Mboi derived for from the two 
methods peaks around ~ 0.1 mag which is not that signif- 
icant. As mentioned in the previous section, we have used 
the average of DENIS & 2MASS A" s -band and MSX D-band 
& ISOGAL 15 /im magnitudes when available. This helps in 
reducing the uncertainties in Mboi and M determinations by 
about a factor \/2. At the same time, it is to be noted that 
large amplitude variable stars such as Miras do show vari- 
ations in the K B -ha.nA upto ~ 2.5 mag (Groenewegen 2006; 
Glass et al. 2001). However, the amplitude of intensity vari- 
ation depends on the period of these sources and is typically 
in the range 0.5 - 1.0 mag. Hence, using only single-epoch 
measurements give us larger errors in the bolometric magni- 
tude determination. The uncertainties due to this variability 
can only be reduced by obtaining time-consuming light curve 
measurements and this is beyond the scope of this present 
work. It is also important to note that there are 239 ISO- 
GAL 15 fim sources which have not been detected in the 
MSX D-band (15 /jm) in our sample. These sources have 
also been used in our sample to determine the Mboi and M. 
The M bol and luminosity (L(L Q ) = i -( M boi-4.75)/2. 5 ) val _ 
ues of each source are given in Table 33, a sample of which is 
given in this paper. Figure [6] shows the histograms of bolo- 
metric magnitudes and luminosities of the sources in MSX 
bulge fields. The approximate bulge RGB tip is predicted 
to be at Ko ~ 8.0 mag (Frogel et al. 1999) which translates 
to M bol ~ -3.5; L « 2000 L . This implies that the distri- 
butions are clearly incomplete below the luminosity of the 
bulge RGB tip. 



6 MASS-LOSS RATE (M) IN THE BULGE 

We have used the dust radiative transfer models for oxygen- 
rich AGB stars from Groenewegen (2006) to derive the mass- 
loss rates of the sources in the MSX bulge fields. These 
models are simulated for values of L = 3000 Lq , distance = 
8.5 kpc, expansion velocity = 10 km s _1 , dust-to-gas ratio 
= 0.005, and no interstellar reddening. Our derivations are 
based on the model for an oxygen-rich AGB star with T e ff of 
2500 K and 100% silicate composition is used. We have used 
the empirical relation between M and (Jf s -[LW3])o to deter- 
mine the mass-loss rate of each star. The luminosity distribu- 
tion of the sources in our MSX bulge fields peaks at ~ 8000 
Lq (see Fig. [6]). Using the scaling law given in Groenewegen 
(2006), we have scaled the model curve for this peak lumi- 
nosity. It should be noted that the empirical values tabulated 
by Groenewegen (2006) are for M < 2 x 10~ 5 M Q yr" 1 . For 
the high mass-loss end (see Fig. 8) we have extrapolated the 
above empirical relation used by us. It should also be noted 



4 Catalogue of MSX sources from the bulge fields with M > 
3 X 10~ 7 Mq yr — 1 is available in the electronic form via VizicR 
Service at the CDS. 
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Figure 6. Histograms showing bolometric magnitudes and lu- 
minosities of the sources in MSX bulge fields which have good 
quality (quality flags 3 and 4) D-band magnitudes and i4" s -band 
photometry. The dotted vertical lines indicate the approximate 
bulge RGB tip (M ho] ~ -3.5; L ~ 2000L©). 



here that the mass-loss rate derivations rely on the accu- 
rate estimates of the extinction and the corrections applied 
to the K s magnitudes. On the other hand, the colour such 
as ([A]-[D])o, which depends little on extinction, allows a 
direct measure of the mass-loss rates. However, we do not 
use this colour for mass-loss estimates primarily because it 
is less sensitive especially for the large and small mass-loss 
rates. The model based on the synthetic colours of AGB 
stars by Jeong et al. (2002) has also been used in litera- 
ture for derivation of mass- loss rate (Ojha et al. 2003). We 
show a comparison of the Groenewegen's model with that of 
Jeong et al. (2002) which is significantly different. In Fig. [7] 
we present the mass-loss rates derived as a function of (K s - 
[LW3])o, using the empirical model of Groenewegen (2006) 
and the empirical relation (M vs. (K s -15)o) of Jeong et al. 
(2002). Here, in place of these colours, we use (K s -[D])q. 
This is justified because the central wavelength of MSX-D 
band is 14.65 (im which is close to 15 /im and falls within 
the LW3/ISO filter (12 - 18/Ltm). It is seen that for Vega 
the difference in magnitudes ([D] — [15]) is 0.04 mag which 
is negligible. It is to be noticed that in the whole range of 
mass-loss we are interested in, for M > 10 -7 M© yr~\ the 
curve based on Groenewegen's model gives lower values of 
mass-loss by factors > 3 as compared to the scale of Jeong et 
al. (2002) for similar colour indices. In this paper, we present 
the mass-loss rates derived from the models of Groenewegen 
(2006). The two mass-loss rate scales are also presented later 
in the discussions of the nature of the MSX bulge sources 
(see Figs. IIT1 and Il2j). 

Fig. [8] shows the if s -[D]o colour distribution of MSX 
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Table 3. Sample catalogue of MSX sources from the bulge fields with M > 3 X 10 7 Mq yr 1 . Data for each source are displayed in two 
lines, with the MSX standard name (e.g. MSX6C-G359. 4989-03. 301) and the IRAS name if any. In order to make the easy comparison 
with available ISOGAL 7 and 15 fim magnitudes, [7] & [15], and IRAS flux densities (in Jy) at 12 & 25 fjxa, S12 k. S25, MSX intensities 
are given in magnitudes for bands A (~ 8 fira) and D (~ 15 fJ.ro), and in flux densities (in Jy) for bands C (~ 12 fim) and E (~ 21 fim). 
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Figure 7. The figure shows the mass-loss rates as a function of 
(fC s -[D])o colour (see text) for the Groenewegen (dotted line) and 
Jeong (solid line) models. 

sources in the "inner" and "outer" zones in the bulge fields. 
There are certain issues regarding the mass-loss estimates 
which need to be mentioned. It should be noted that using 
the empirical model of Groenewegen (2006) includes uncer- 
tainties in the determination of the mass-loss rates arising 
from the following reasons: (1) The relation was derived for 
oxygen-rich AGB stars in the solar neighbourhood. It has 
also been argued that metallicity affects the dust-to-gas ra- 
tio and the outflow velocity from evolved stars (Habing et al. 
1994). This is directly related to the mass-loss and, there- 
fore, the colour-mass-loss relation could possibly differ in 
different environments such as between the Galactic bulge 
and the Magellanic Clouds. (2) It is important to emphasize 
that the AGB stars are variable in nature. The average K 



amplitude of sources associated with known LPVs (Glass et 
al. 2001) is ~ 1.0 mag, which amounts to a factor of ~ 2 
- 5 uncertainty in the determination of M. (3) The model 
(Groenewegen 2006) used is valid for a limited range of lu- 
minosity and expansion velocity and are simulated for AGB 
stars. Groenewegen (2006) have also mentioned about the 
caveat of post-AGB model simulation. Here the model for 
post-AGB sources is calculated under the assumption that 
the effective temperature and luminosity do not change over 
the time for the dust shell to drift away. So using the colour- 
mass-loss relation for sources having these parameters be- 
yond the range (i.e sources with high luminosity and expan- 
sion velocity like PNe and post-AGB or sources with lower 
luminosity like T-Tauri stars) will have a large uncertainty 
in their mass-loss estimation. The classical AGB luminosity 
limit is LogL(L Q ) ~ 4.74 (Marigo et al. 1998; Zijlstra et 
al. 1996). Hence, in our MSX bulge fields, sources having 
luminosity values beyond this cut-off limit will have larger 
uncertainties in the mass-loss estimation. These sources and 
their contribution to the integrated mass-loss rate will be 
discussed in detail later in this section. 

The mass-loss rates for the objects in the MSX bulge 
fields range from 10~ 7 to 10~ 4 Mq yr -1 . As is obvious from 
the comparison of the two histograms in Fig. |SJ the number 
of sources is incomplete for M < 3 x 10~ 7 MQyr _1 (Groe- 
newegen's model) in the "inner" zone, and, more severely, 
for M < 2 x 10~ 6 M yr" 1 in the "outer" zone. The source 
excess seen particularly in the "inner" zone for jf s -[D] ~ 
9 - 10, is most likely due to spurious association of the 
sample-A and sample-U sources (see §2), which makes the 
estimated mass-loss rates appear smaller than their actual 
values. Numerical values of the mass-loss rates of individ- 
ual sources are given in Table 3. The number distribution 
of mass-loss rates in the "inner" and "outer" bulge zones 
as a function of M are displayed in Fig. [9] In Fig. 1101 we 
plot the corresponding values of the average total mass-loss 
rate per square degree and per 0.5 bin of logM for the "in- 
ner" zone. We do not attempt to build the same plot for the 
"outer" zone because all the mass-loss bins in this zone, ex- 
cept M ~ 0.3 — 3 x 10 Mq yr" 1 , would be very uncertain. 
This is due to the incompleteness of the lower mass-loss bins 
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Figure 9. Distribution of mass-loss rates (M) of MSX sources in the "inner" (left) and "outer" (right) bulge zones. The mass-loss rates 
are inferred from the Grocnewegen's model using the (K B -[D])o colour. 
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Figure 8. i^ s -[D] colour distribution of MSX sources. The mass- 
loss rate scales displayed at the top of the figure are from the 
models of Grocnewegcn (upper panel) and Jcong (lower panel). 
The full line histogram represents the stars of the "inner" zone 
< 0.5°) and the dashed line denotes stars of the "outer" zone 
(0.5° < \l\ < 3.0°) in the bulge (Fig.[2]and Table 1). 
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Figure 10. Average total mass-loss rate per square degree and 
per 0.5 bin of log (M) as a function of M for MSX sources in the 
"inner" bulge zone. 



and the uncertainty in K S -[D] for M > 3 x 10" 5 M© yr" 1 . 
Numerical values of the integrated mass-loss rate in the "in- 
ner" zone are displayed in Table 4, where we have limited 
the integration to mass-loss rates, M > 3 x 1CP 7 M© yr^ 1 
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for which the data are reasonably complete. The first row 
shows the value of mass-loss rate per square degree averaged 
over all the fields in the "inner" zone. In the next five rows, 
we give the integrated mass-loss in the five latitude bins 
which cover the bulge zone. The table also lists the corre- 
sponding integrated mass-loss rate per unit stellar mass (in 
yr" 1 ). The integrated mass-loss rate in the entire "inner" 
zone for the bulge fields is 1.96 x 10 -4 Mq yr -1 deg -2 and 
the corresponding integrated mass-loss rate per unit stellar 
mass is 0.48 x 10 -11 yr -1 (see Table 4). For the rest of the 
discussion that follows, it should be kept in mind that the 
integrated mass-loss rate is only for the "inner" zone and 
the contribution from the "outer" zone is not included. 

Here, we would like to discuss about the sources with 
luminosities beyond the classical AGB limit of LogL(L,Q) > 
4.74. There are 18 such sources in the "inner" zone for which 
we derive the integrated mass-loss rate. Out of these only 
three sources have mass-loss rates M > 3 x 10 -7 Mq yr -1 . 
The uncertainty in the mass-loss estimation of these three 
sources will not affect the integrated mass-loss rate as they 
contribute less than 4% to the total mass-loss. We have 
also examined the nature of all these 18 sources in vari- 
ous CM and CC diagrams. Except for one source which has 
a faint As-band magnitude (12.3 mag), rest of sources are 
very bright in K s (< 5.6 mags) with 4 of them saturated 
in 2MASS (K a < 3.5 mag). All of these sources are also 
bright in the MSX D-band (< 2.8 mag). Referring to Figs. [TT] 
and 1131 we see that the majority of these sources have blue 
colours ((AV[D]) ~ 2; ([A]-[D]) ~ 1). This implies that 
they are mostly foreground sources which have possibly not 
been removed using our foreground source rejection proce- 
dure (see §3). The foreground (or early- type) nature of these 
sources were further verified in the (I — J)/(K S — [D]) CC 
diagram based on the discussion presented in Schultheis et 
al. (2002; cf Fig. 2). Among these 18 sources, there are two 
sources which show large colour excess in the diagrams dis- 
cussed above. These are likely to be YSOs. Additional spec- 
troscopic and photometric observations are required to un- 
derstand the nature of these luminous sources. However, the 
uncertainties involved in the mass-loss rates of these sources 
do not influence the estimation of the integrated mass-loss 
rate as their contribution is negligible as is already men- 
tioned. 

We have also calculated the integrated mass-loss rate 
for the "intermediate" (|/| < 3.0°, 1° < \b\ < 2°) and 
the "outer" < 3.0°, 2° < |6| < 5°) bulge fields. For 
M > 3 x 10 M©yr , the integrated mass-loss rates are 
3.87 x 10- 4 M Q yr- 1 deg- 2 and 1.32 x 10~ 4 M Q yr" 1 deg" 2 
in the two bulge fields, respectively. In comparison, Ojha 
et al. (2003) derive values of 3.7 x 10~ 4 M Q yr" 1 deg" 2 and 
0.4 — 1.0 x 10" 11 yr" 1 for the integrated mass-loss rate and 
mass-loss rate per unit stellar mass, respectively, using ISO- 
CAM 7 and 15 /im observations. It should be noted here that 
the Galactic bulge fields presented in Ojha et al. (2003) cover 
-1.5° < I < +1.6°; -3.8° < b < -1.0°, b = +1.0° with a 
total area of ~ 0.29 deg 2 whereas the area covered in this 
present work is much larger and hence offers better statistics 
than the sample of Ojha et al. (2003). The sensitivities of 



Table 4. Integrated mass-loss (in Mq yr - 1 deg -2 ) and integrated 
mass-loss rate per unit stellar mass (in yr" 1 ) in the "inner" bulge 
zone with M > 3 X 10 _7 MQyr~ 1 . The mass-loss rates are cal- 
culated using the model of Groenewegen (2006). 
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the MSX and 2MASS data also allow us to probe the high 
mass-loss end with a better statistical sample. 

Le Bertre et al. (2001; 2003) have studied the Galac- 
tic mass-losing AGB stars and concluded that the replen- 
ishment to the ISM is dominated ( > 50%) by AGB stars 
with mass- loss rates > 10 -6 Mq yr -1 . These sources consti- 
tute ~10% of their sample and noticeably there are no AGB 
stars in their sample with mass-loss rates > 10 _4 MQyr _1 
though sources with mass-loss > 10 _4 MQyr _1 are known 
to exist (Habing 1996). The sensitivities of the data sets 
used by us probe these rare high mass-loss stars. In our MSX 
bulge fields, we have a total of 42 sources displaying mass- 
loss > 1.0 x 10~ 5 Mq yr -1 which significantly contribute to 
the mass returned to the ISM. In comparison, the stellar 
population studies in the solar neighbourhood by Jura & 
Kleinmann (1989) identify 63 sources with mass- loss rates 
> 10" 6 Mq yr" 1 out of which only 21 sources in their sam- 
ple have mass-loss rates > 1.0 x 10" 5 Mq yr" 1 . From their 
sample it is seen that oxygen- and carbon-rich AGB stars 
have almost equal contribution to the replenishment of the 
ISM which amounts to 3 - 6 x 10" 4 M kpc -2 yr" 1 . From 
the MSX bulge fields we estimate the total mass replenish- 
ment to the ISM to be ~ 9 x 10" :! M kpc" 2 yr -1 , which 
is not surprising given the large number of high mass-loss 
AGB stars (M > 1.0 x lO^Meyr" 1 ) detected in our sam- 
ple. However, it should also be kept in mind that the sample 
of Jura & Kleinmann (1989) is incomplete at both the lower 
and the high mass-loss ends. Since the low mass-loss end 
is incomplete in our sample, we refrain from making any 
comparison with other studies. 



7 NATURE OF THE MSX SOURCES 

In this section, we discuss the various colour-magnitude and 
colour-colour diagrams. This enables us to study the nature 
of the sources in the MSX bulge fields. As has been men- 
tioned earlier, we present only sources with good quality 
MSX and 2MASS data. 
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7.1 Colour-magnitude diagrams 

In Figs. 1111 1121 1131 and 1141 we present the dereddened 
(JC s -[D])o/[D] , (K s -[V]) /K s0 , ([A]-[D])o/[D] and (K s - 
[A])o/[A]o colour- magnitude diagrams of MSX sources in 
the two bulge zones, respectively. 

Figures [TT] and [TJ] show the two colour- magnitude dia- 
grams involving the (-?Cs-[D])o colour which has been used 
in §6 for the derivation of the mass-loss rates of the sources. 
The mass-loss rate scales derived from Jeong et al. (2002) 
and Groenewegen (2006) are displayed on the top panel of 
these two figures. It is interesting to see that the sources from 
sample-A (which are plotted in red colour) are mostly high 
mass-loss stars with (AT s -[D])o > 6.5, corresponding to a 
mass-loss rate > 5 x 10~ 6 Mq yr _1 based on the Groenewe- 
gen's scale. The sample-U sources (shown in blue) extend 
the high mass-loss sequence beyond the sample-A sources. 
It should be kept in mind that the straight line sequence seen 
for the sample-U sources is partly due to the assigned lower 
limit of 13.7 mag on the 7f s -band magnitude. It is worth 
recalling here that these two samples constitute sources 
brighter than 4th mag in MSX D-band with no 2MASS 
counterpart with K s < 11.0 mag within 4" search radius. 
Sample-A includes sources having a 2MASS counterpart 
within an extended 5" search radius without any _K" B -band 
magnitude cut and sample-U comprises of sources having no 
2MASS counterpart within this radius. We have assigned K s 
= 13.7 mag for all the sample-U sources. Another distinct 
feature which is seen in the "inner" zone is branching of the 
sources into two sequences with an appreciable gap around 
its-[D])o > 5.5. It is also interesting to note that majority 
of the ISOGAL sources seen in the "inner" bulge zone, i.e. 
lower luminosity sources, occupy the low mass-loss end of the 
plots. Our prime focus in this study is to identify the nature 
of high mass-loss stars in the bulge fields. We therefore se- 
lected sources with large mass-loss rates (> 3x 10~ 5 MQyr -1 
as per the Groenewegen's scale) and surveyed the VizieR 
and SIMBAD database for identification with known ob- 
jects. Since majority of these high mass-loss sources are 
known IRAS sources, we used a search radius of 10" by 
taking into consideration the typical IRAS PSC error ellipse 
of ~ 10" x 20". We have grouped these sources identified in 
VizieR/Simbad into four major classes - Planetary nebulae 
& post AGB stars, maser & OH/IR sources, peculiar sources 
(which also include variable stars) and IRAS sources which 
only have IRAS names but no other classifications. We have 
marked these four classes in the colour-magnitude diagrams 
using different symbols (see caption of Fig. Ill[) . 

In Fig.[l3l we have plotted the ([A]-[D]) /[D] colour- 
magnitude diagram. There is a gap seen around ([A]-[D])o 
~ 2.5. Among the sources to the right of this gap (the 
redder sources) the majority are known post AGB and 
planetary nebulae. In this figure it is also clearly evident 
that sources from sample-A are redder than the rest. Fig- 
ure !14l shows (_ftr a -[A])o/[A]o colour-magnitude diagram. The 
(i^s-[A])o/[A]o diagram seems to be the best criterion for 
the detection of large amplitude LPVs (Glass et al. 1999, 
Schultheis et al. 2000). We see two interesting sequences in 
this colour-magnitude diagram. The sources seem to branch 
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Figure 11. [D]o/(Ks-[D])o magnitude-colour diagram for the 
sources in the MSX bulge fields. The upper panel shows the 
"outer" zone (0.5° < |i| < 3.0°) and the lower panel shows 
the "inner" zone (|/| < 0.5°) of the bulge. Only good quality 
data (MSX flags 3 & 4; 2MASS 'rd-flg' 1-3) are presented. 
The colour coding which is followed is red for sample-A sources, 
blue for sample-U sources and black for rest of the bulge fields. 
The mass-loss rate scales displayed on top are from Groenewegen 
(2006) and Jeong et al. (2002) (lower) models. The two curves 
are from the models of Groenewegen (2006), where the red curve 
represents the model for an oxygen-rich AGB star with T e ff = 
2500 K and a 100% silicate composition. The blue curve is for 
a carbon-rich AGB star with T e g = 2650 K with 100% amor- 
phous carbon (AMC). The model curves are scaled to the peak 
luminosity of 8000 Lq of the sources in the bulge fields. In case 
of the Groenewegen models, the magnitude in the ISOCAM LW3 
(15 (im) filter has been used in place of the MSX D-Band. The 
ISOGAL counterparts are marked as magenta open squares for 
the entire bulge fields including the samplc-A sources. The known 
sources from the VizieR and SIMBAD database are displayed 
with the following symbols while maintaining the colour coding 
for the sample-A, sample-U and the rest of the bulge fields sources 
: Planetary Nebulae and Post AGB stars - Solid squares; Maser 
and OH/IR sources - solid triangles; Peculiar sources including 
variable stars - ffi; IRAS sources with no other classifications - 
crosses. The two carbon stars in our sample are represented by 
0. Sources saturated in the 2MASS fC s -band (< 3.5) are shown 
as stars. 



out into two for (if s -[A])o > 2, with an appreciable gap 
between the two. The upper sequence primarily consists of 
sources with J — K excess and the lower sequence mostly 
comprises of sources which do not have any D-band associ- 
ation. Interestingly, the majority of the sources in the lower 
sequence do not have either C- or E-band association and 
are fainter in K s (> 9 mag). The sources from sample-A 
seem to extend the sequence with increasing slope of the 
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Figure 12. K a o/(K B -[D])o magnitude-colour diagram of MSX 
sources with 2MASS counterparts. The symbols are same as 
shown in Fig. 1111 




Figure 13. [D]o/([A]-[D])o magnitude-colour diagram of good 
quality MSX sources. The IRAC 8 /an magnitude is used in place 
of MSX A-band for the Groenewegen models. The symbols are 
same as shown in Fig. 1111 



curve with increasing (Jf s -[A])o colour and the majority of 
these are again known post AGB and planetary nebulae. 
The sample-U sources seem to form a parallel sequence and 
occupy the red end. As has been mentioned earlier, it must 
be noted that this sequence is also a result of the assigned 
lower limit of 13.7 mag on the _K" s -band magnitude of the 
sources in this sample. In both the above diagrams the ISO- 
GAL sources are towards the blue end. 



7.2 Colour-colour diagrams 

The colour-colour diagrams are also useful for discriminat- 
ing between different classes of objects and to study their 
nature. Ortiz et al. (2005), have studied the evolution of 
carbon- and oxygen-rich AGB stars, post-AGB stars and 
planetary nebulae using mainly data from MSX. In compar- 
ison to their Fig. 3, we have plotted the (D — E)/(A — D) 
colour-colour diagram in Fig. 1151 We have plotted the red- 
dened colours to facilitate comparison with the plot of Ortiz 
et al. (2005). As noted by Ortiz et al. (2005), we also find a 
distinct gap in the plot which is shown as a solid line based 
on the visual inspection of our data, in addition to the line 
(dotted one) from Ortiz et al. (2005). The solid line based on 
our estimate allows to discriminate few additional border- 
line objects. In this figure, we have plotted the known high 
mass-loss sources, which have already been identified. Apart 
from these we have also searched the VizieR and SIMBAD 
database for known counterparts of all the additional sources 
lying to the right of the gap, which were missed out in the 
earlier selection based on Fig. 1111 which included only good 
quality its-band sources. Majority of the known IRAS and 
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Figure 14. [A]o/(A" s -[A])o magnitude-colour diagram of MSX 
sources (quality flags of 3 & 4) with 2MASS counterparts. The 
upper panel shows the "outer" zone (0.5° < \l\ < 3.0°) and the 
lower panel shows the "inner" zone < 0.5°) of the bulge. The 
symbols are same as in Fig. 1111 The cyan circles represent the 
sources without D associations. 
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Figure 15. (D — E) / (A — D) colour-colour diagram for the two 
bulge zones are presented (open black circles). The known ob- 
jects from the VizieR and SIMBAD database are displayed. The 
symbols are same as shown in Fig. 1111 The dotted line from Or- 
tiz et al. (2005) shows the position of the gap between. We have 
added the solid line which allows one to identify a few additional 
border-line objects 



maser sources lie to the left and lower end of the plot. Plan- 
etary nebulae populate mostly the region beyond the gap. 
Ortiz et al. (2005) suggest that the gap seen is due to the 
rapid evolution of stars as they cross the dotted line which 
can be considered as a transition boundary. This is in good 
agreement with the different population of sources found on 
either side of the gap in Fig. 1151 

Using MSX and 2MASS mid- and near-infrared colour- 
colour diagrams, Lumsden et al. (2002), Sevenster (2002) 
and Messineo et al. (2004) have studied different classes of 
evolved stars with circumstellar envelopes. In Fig. 1161 we 
present the different colour-colour diagrams of the MSX 
bulge fields for good quality flags in all the MSX bands 
and also in the iis-band. The plots are shown in the units 
of the reddened flux ratio rather than the colour index for 
easy comparison with the plots of Lumsden et al. (2002) . In 
Fig. 116b . we plot the F21 /Fg vs. Fg /Fk flux ratios. Accord- 
ing to Lumsden et al. (2002), the oxygen- and carbon-rich 
populations separate out clearly in this colour-colour dia- 
gram (see Figs. 10 & 11 of Lumsden et al. 2002). However, 
we do not see any such separation here in our MSX bulge 
fields. The OH/IR stars and the planetary nebulae mostly 
have F21/F8 flux ratios > 1 consistent with the trend seen 
in Fig. 10 of Lumsden et al. (2002). The mid-infrared colours 
of planetary nebulae are similar to HII regions but mostly 
planetary nebulae are bluer because they are not embed- 
ded in molecular clouds. Figure [T&b shows the Fi2/Fg vs. 
F21/F14 colour-colour diagram. This diagram is useful to 



distinguish between the AGB and the post-AGB phases and 
to locate post-AGB stars which have redder colours due to 
their colder envelopes (Sevenster 2002; Messineo et al. 2004). 
The transition from a blue (< 2.38) to red (> 2.38) Fi 2 /F 8 
flux ratio may correspond to a transition off the AGB to 
proto-planetary nebulae and from a blue (< 1.90) to red (> 
1.90) F21/F14 flux ratio indicates a later evolutionary tran- 
sition, when mass-loss starts to drop down several order of 
magnitudes and there is the onset of the fast wind (Messineo 
et al. 2004). Most of the MSX sources show Fi 2 /F 8 < 2.38 
and F21/F14 < 1.90 as expected for AGB stars. In Fig. [Tot, 
we plot the F21/F14 vs. F14/F8 colour-colour diagram. We 
do not see any trend separating the oxygen-rich population 
from the carbon-rich one as pointed out by Lumsden et al. 
(2002). The majority of the sources are OH/IR stars and 
IRAS sources which occupy the central region of the plot. 
The planetary nebulae are mostly above F2i/Fi4>l and 
Fi4/Fg>l. We see a gap around Fi4/Fg ~ 3. This could 
possibly suggest the transition from the AGB to the more 
evolved phase. In Fig. I16H. we plot the F21/F8 vs. F14/F12 
flux ratios. The MSX sources with F 2 i/Fg > 5.96 (see also 
Fig. [T6k ) probably represent post-AGB stars (Messineo et 
al. 2004). 

7.3 Population of carbon stars in the bulge fields 

We have investigated the nature of the red sources (with 
J — K s > 3) in the various diagrams presented. They oc- 
cupy the high mass-loss end of Figs. llll and [T2l as expected. 
From the loci of the Groenewegen's models, it is difficult 
to comment on whether these sources are carbon-rich or 
oxygen-rich AGB stars. The population of carbon-rich stars 
if present in the bulge would be crucial in our understand- 
ing of the bulge formation, star formation history in the 
bulge and the nature, formation and evolution of this rare 
population itself. Except for the peculiar 34 carbon stars de- 
tected by Azzopardi et al. (1991), there has been no other 
carbon star detection in the bulge. The bulge membership 
and nature of these stars have been a matter of much de- 
bate (Ng 1998; Whitelock 1993). In Figs. [11] and [12] we 
see sources falling on the locus of the carbon-rich AGB star 
model. But based on just the colour and magnitude, it is im- 
possible to confirm or reject the possibility of finding carbon 
stars in the bulge. Apart from observations focussed on the 
Galactic bulge, there are several other studies based on the 
near- and mid-infrared colour-colour and colour-magnitude 
diagrams to identify carbon star populations in different 
Galactic environments (van Loon et al. (1998); Buchanan 
et al. (2006); Thorndike et al. (2007); Groenewegen et al. 
(2007)). van Loon et al. (1998) have presented the (H - K) 
vs (K — [12]) colour-colour diagram to distinguish between 
carbon and oxygen rich AGB stars. The Figs. 3a and b of 
van Loon et al. (1998) show the Galactic sample of AGB 
stars from Guglielmo et al. (1993). In their figures, the 
oxygen and carbon rich sequences clearly separate out. In 
Fig. 1171 we have plotted the similar colour-colour diagram 
((H— K s )/(K a — C)) of the sources in our bulge fields. In this 
plot we have plotted only sources from sample-A and rest of 
the bulge fields. Sources from sample-U, which do not have 
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Figure 16. The near- and mid-infrared colour-colour diagrams for the two bulge zones (open black circles). The rectangular box shown 
in (b) is based on the criteria discussed in Messineo et al. (2004). The symbols are same as shown in Fig. 1111 



_ff-band magnitudes, are excluded. To compare with the fig- 
ures of van Loon et al. (1998), we have plotted here the 
reddened colours. The scatter in our plot makes it difficult 
to differentiate the two sequences. In their study of lumi- 
nous sources in the LMC, Thorndike et al. (2007) have used 
the ((H — K)/(K — [A])) colour-colour diagram to classify 
different stellar populations. In Fig. 1181 we have plotted a 
similar reddened ((H — K a ) / (K s — A)) colour-colour diagram 
for the sources in the MSX bulge fields. Here also we see the 
branching off sources into two sequences similar to that seen 
in Fig. 1141 The oxygen- and carbon- rich AGB star models of 
Groenewegen (2006) trace the sources in the upper sequence 
better. However, given the scatter in the plot, it is difficult 
to identify any carbon-rich AGB star sequence as shown in 
Thorndike et al. (2007). We also derive similar conclusions 
from the investigation of the ((J — K a )/(K a — [A])) colour- 
colour diagram (not presented in this paper) based on the 
plots presented in Buchanan et al. (2006) and Groenewegen 
et al. (2007). However, it is interesting to note that there 
are two known carbon stars (IRAS 17534-3030 and IRAS 
17547-3249) in our sample, one each in the "inner" and the 
"outer" zone. These sources classified as carbon stars (Groe- 



newegen et al. 2002; Volk et al. 2000; Guandalini et al. 2006) 
do not seem to fall on the carbon-rich model curve shown in 
Figs. HU HU [13] and [TJ] One of the two carbon stars (IRAS 
17534-3030), which lies in the "inner" zone of sample-A, has 
good quality 2MASS H, K s , MSX A- and C-band data and 
is shown in Fig. [I7]and Fig. 1181 It is seen to lie far away from 
the carbon line. Here, we would like to recall that for sample- 
A sources, the chance associations are large (~ 0.6; see §2) 
and hence the possibility of a spurious 2MASS association. 
However, for this particular source which is from sample-A, 
we have checked the near- and mid-infrared colours in detail. 
Comparing the colours with the plots presented in Lumsden 
et al. (2002), we see that the flux ratios agree well with the 
carbon star population except for the F2i/Fg and F21/F12 
ratios which are marginally on the redder side. It is worth 
mentioning here that this carbon star IRAS 17534-3030 is 
a well known extreme carbon star which has been recently 
studied by Pitman et al. (2006) using ISO-SWS spectra. 
They predict the presence of nitride dust in the star. Apart 
from this, only near-infrared colour-colour diagrams have 
also been used to identify carbon stars. Based on the study 
of Kerschbaum & Hron (1994), we have also checked the 
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Figure 17. (H — K B )/(K B -[C]) colour-colour diagram for the "in- 
ner" (solid circles) and "outer" (open circles) zones are presented. 
The sources from samplc-A are shown in red. The carbon-rich (red 
line) and oxygen-rich (blue curve) sequences from van Loon et al. 
(1998) are also plotted. One of the two carbon stars identified in 
our sample is marked with 0. 

(J — H) vs (H — Ks) colour-colour diagram (not presented 
in this paper). Here also, we are unable to identify any car- 
bon star sequence except for the identified carbon star IRAS 
17534-3030 whose (J — H) / (H — Ks) colours are consistent 
with the locus of carbon stars. Lebzelter et al. (2002) and 
Cioni & Habing (2003), have shown the use of the (/ — J) 
vs (J — K) colour-colour diagram to discriminate between 
the oxygen- and carbon-rich populations for the LMC. Since 
most of the sources (including the two known carbon stars) 
in our sample do not have / counterparts, we have not pre- 
sented this colour-colour diagram in this paper. From the 
figure of Ortiz et al. (2005) , the carbon stars are seen to pop- 
ulate the lower left corner (bluer end - [14.7]-[21.3] colour 
ranging from ~ 0.2 to 0.8 and [8.3]-[14.7] colour ranging from 
~ to 1.5) of the diagram. One of the two carbon stars in 
our sample falls in this region (see Fig. I15|l . Other than this, 
there is no such distinct population seen among the sources 
belonging to our bulge fields. We infer that based on only 
the colour-colour and colour magnitude diagrams we cannot 
conclusively comment on the presence or absence of carbon 
stars in the bulge. 



8 CONCLUSION 

In this paper, we have presented in detail the study of the 
AGB population detected by MSX in the "intermediate" 
(|Z| < 3.0°, 1° < \b\ < 2°) and "outer" < 3.0°, 2° < 
|6| < 5°) Galactic bulge fields covering a large area of 48 
deg 2 . We have shown that the MSX data in conjunction 
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Figure 18. (H — K S )/(K S -[A]) colour-colour diagram for the "in- 
ner" (solid circles) and "outer" (open circles) zones are presented. 
The symbols are same as shown in Fig. 1111 



with the 2MASS database can be potentially used to detect 
AGB stars well above the RGB tip with the determination 
of their luminosities (provided they belong to the bulge) and 
mass-loss rates. The sensitivities of the two surveys enable us 
to sample the high mass-loss end, characterized by excess in 
A' s -[D])o colour due to circumstellar dust emission, with bet- 
ter statistics. We have derived the mass-loss rates of all the 
MSX sources which range from 10~ 7 to 10~ 4 Moyr _1 and 
estimated the integrated mass-loss rate in the MSX bulge 
fields. Taking the completeness into account, we have limited 
the integration to mass-loss rates, M > 3x 10 -7 Mq yr _1 for 
the "inner" zone < 0.5°) only. There is a factor of ~ 3 
increase in the estimated integrated mass-loss rate in the 
"intermediate" (3.87 x 10 -4 Mq deg -2 yr _1 ) as compared 
to the "outer" (1.32 x 10~ 4 M Q deg -2 yr _1 ) bulge regions. 
The average integrated mass-loss rate is estimated to be 
1.96x 10 -4 Mq deg -2 yr _1 and the corresponding integrated 
mass-loss rate per unit stellar mass is 0.48 x 10 _ yr . 
Apart from the mass-loss derivations, we have used the var- 
ious colour-magnitude and colour-colour diagrams to discuss 
in detail the nature of the MSX sources in the bulge fields 
and identify the location of planetary nebulae, post-AGB 
stars, OH/IR sources etc. in these diagrams. 

Studies based on proposed SpitzerlRAC and MIPS pho- 
tometric surveys of the innermost Galaxy including the most 
obscured and crowded region would in future enable us to 
unravel in detail the late stages of stellar evolution partic- 
ularly of the infrared-luminous AGB stars, the mass-loss of 
which is crucial in understanding the later stellar fate and 
hence the dust in the universe. Future infrared mission (e.g. 
WISE ; Duval et al. 2004) hold the key to such studies. 
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